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Abstract 



o 

CNj ' If gravity is an emergent phenomenon, as suggested by several recent results, then the structure 

, of the action principle for gravity should encode this fact. With this motivation we study several fea- 

tures of the Einstein-Hilbert action and establish direct connections with horizon thermodynamics. 
We begin by introducing the concept of holographically conjugate variables (HCVs) in terms of which 
the surface term in the action has a specific relationship with the bulk term. In addition to g ab and 
its conjugate momentum y/—gM cab , this procedure allows us to (re)discover and motivate strongly 
the use of f ab = y/—gg ab and its conjugate momentum N^ b . The gravitational action can then be 
' . interpreted as a momentum space action for these variables. We also show that many expressions 

in classical gravity simplify considerably in this approach. For example, the field equations can be 
written in the form d c f ab = &H g /dN^ b ,d c N^ b = —dH g / 'df ab (analogous to Hamilton's equations) 
j^Q' for a suitable Hamiltonian ~H g , if we use these variables. More importantly, the variation of surface 

term, evaluated on any null surface which acts a local Rindler horizon can be given a direct ther- 
modynamic interpretation. The term involving the variation of the dynamical variable leads to T5S 
while the term involving the variation of the conjugate momentum leads to S8T. We have found this 
correspondence only for the choice of variables (g a b, \J—gM cab ) or (f ab , N^ b ). We use this result to 
provide a direct thermodynamical interpretation of the boundary condition in the action principle, 
when it is formulated in a spacetime region bounded by the null surfaces. We analyse these features 
from several different perspectives and provide a detailed description, which offers insights about the 
nature of classical gravity and emergent paradigm. 



1 Introduction 



— * 

The curious connection between gravity and thermodynamics first came to light with the work of Beken- 
stein [1-3], which ascribed to a black hole an entropy proportional to the surface area of its horizon. 
Soon, it was discovered that black hole horizons possess temperature as well [4,5]. In the four decades 
since then, the intriguing connection between gravity and thermodynamics has been steadily becoming 
stronger (see e.g., [6,7] ). 

One paradigm which has emerged from this connection considers the dynamics of gravity as not of 
fundamental nature, but emergent from the dynamics of a more fundamental theory, just as the thermo- 
dynamics of a material system emerges out of the basic dynamics of its molecules. (For a recent review, 
see [8,9]; for other work similar in spirit, see e.g., [10-14]). This paradigm has obtained support from 
the following results: 

a) Gravitational field equations in a wide class of theories - more general than Einstein gravity - lend 
themselves to a thermodynamical interpretation [15-17]; 

b) Gravitational equations of motion can be obtained from thermodynamical extremum principles [18, 
19]; 
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c) It has been possible to obtain the density of microscopic degrees of freedom through equipartition 
arguments [20,21]; 

d) The action functional for gravitation in a class of theories has a thermodynamic interpretation [15, 



e) Einstein's equations reduce to Navier-Stokes equations of fluid dynamics in any spacetime when 
projected on a null surface [25,26]. This generalizes previous results for black hole spacetime [27-29]. 

f ) The euclidean path integral of the gravitational action interpreted as a partition function has provided 
expressions for free energy, energy and entropy [30] in the Lanczos-Lovelock models generalizing 
previously known results [15,31]. 

These results show that gravity could indeed be an emergent phenomenon, similar to elasticity or fluid 
mechanics. 

Conventionally, however, one treats gravity like any other field and its dynamics is obtained from a 
standard action principle. The fact that the theory obtained by such a variational principle possesses 
an emergent description strongly suggests that the action functional itself must encode this information. 
The main purpose of this paper is to explore several aspects of gravitational action principle from the 
emergent gravity perspective and unravel these features to the extent possible. 

In fact, it is well-known that general relativity has some peculiarities which introduces special difficulties 
(not encountered in other field theories, e.g. non-Abelian gauge theories) when we try to derive the 
Einstein equations from an action principle (see e.g., Chapter 6 of [32]). There does exist a scalar action, 
the Einstcin-Hilbert action, which, when added to the matter action, can be varied with respect to the 
metric to obtain the Einstein equations. Based on the usual theoretical prejudice [33], the equations 
of motion which are second order in derivatives of the coordinates, are obtained from an action which 
is quadratic in the first derivatives of the dynamical variables. The generally covariant Lagrangian for 
gravity, however, is forced to be at least second order in derivatives of the metric since any scalar made of 
the metric and its first derivatives will vanish in the local inertial frame and, by virtue of being a scalar, is 
zero in any other frame. The usual choice, R, has a special structure (viz. linearity in second derivatives) 
which allows us to obtain the equations of motion which are only second order in the derivatives of the 
metric, if (and only if) we fix metric and its derivatives on the boundary. This is possible because we 
can separate out the the second derivatives in the Lagrangian into a total derivative which becomes a 
surface term in the action. Its variation does not contribute to the equations of motion if we set the 
variations of the metric and its derivatives to be zero on the boundary. 

The main conceptual difficulty with this program — which makes gravitational action different from 
those in other field theories — is that the we need to fix both the dynamical variable and its derivative 
at the boundary to obtain the equations of motion. There are some issues with this procedure. Suppose 
that the spacetime region we are considering is the region between two spacelike surfaces. Assume 
that all quantities go to zero at spatial infinites. If we now fix the metric and its derivative on the 
earlier time-slice, the Einstein equations would give us the corresponding values on the latter time-slice. 
Thus, we do not really have the freedom to fix arbitrary values for the metric and its derivatives on the 
boundaries. Also, setting our eyes on a quantum theory, we would not want to fix both the dynamical 
variable and its derivative at the same spacelike hypcrsurface. Another option would be to add a term 
to the action, the variation of which will cancel the terms with variation of the first derivative. A well- 
known example is the Gibbons-Hawking- York counterterm [31,34] though it is not unique [35]. But this 
entire approach appears a little contrived and the resulting action, for the Gibbons-Hawking- York case, 
although generally covariant becomes foliation dependent. 

There is, however, another aspect to this issue. If we use g a b as the dynamical variables (with an 
associated canonical momenta), then it turns out that the surface variation involves only the variation 
of the canonical momenta. This is a nontrivial result and arises only because the surface and bulk terms 
of the Hilbert Lagrangian are connected in a peculiar manner [22]. This, in turn, implies that one can 
treat the Hilbert action as a momentum space action and the equations of motion can be obtained by 
fixing the canonical momenta at the boundary. Thus, gravity may be better dealt with in the space of 
the canonical momenta corresponding to the metric rather than in the space of the metric. Further, we 
will demonstrate that this program will work with the components of g a b as variables, but will fail with 
g ab as variables. Although g ab is a failure, we shall find that the components of the corresponding tensor 
density, f ab = \ / —gg ab , would also serve our purpose. 




22-24]; 
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The variables f ab = \J~ r gg ab have been used in classical literature [36-39] but somehow does not seem 
to have attracted sufficient attention in recent years. One of the purposes of this paper is to advertise 
nice features of these variables and some pedagogical results. After a slight digression on how various 
expressions simplify on being expressed in terms of f ab and the corresponding canonical momenta, we 
look more closely at the result that an integral of the surface term in Einstcin-Hilbcrt action over a horizon 
gives the entropy [15,22-24,40,41] of the horizon. The results have been obtained for static metrics in 
which the action will be proportional to the range of time integration r due to the static nature of the 
metric. The integrals of the surface Lagrangian over the horizon will give tTS, where r is the range of the 
time integration, T is the temperature and S is the entropy. Usually we work in Euclidean sector where 
the natural choice for r is the inverse temperature ft, giving tTS = S. Alternatively, one can define the 
surface Hamiltonian [40,41] by H SUI = (dA SUI /dr) = TS and concentrate on just the integration over the 
co-ordinates transverse to the horizon. A variation of this integral would provide us with TSS + S5T. We 
will show that the variation of the surface term in the Einstein-Hilbert action for static metrics allows 
us to identify the TSS term as arising from the variation of the dynamical variable g a b or f ab and the 
SST term arises from the variation of the corresponding canonical momenta in both cases. (In contrast, 
such a nice separation docs not work if, e.g., we use g ab as dynamical variables). 

This result extends to a very general class of null surfaces, which act as horizons for local Rindler observers 
and suggests a simple thermodynamical interpretation of the variational principle. In any region bounded 
by null surfaces, one can introduce local Rindler observers who perceive patches of the null surfaces as 
horizons and attribute temperatures to them. Using this, we can reformulate the boundary condition for 
the gravitational action principle as equivalent to keeping the temperature constant on these null surfaces. 
We believe this offers some insight into the thermodynamic interpretation of action principle. 

The outline of the paper is as follows. In Section 2, we review some well-known results on the structure 
of the Einstein-Hilbert action to set the stage. In Section 3, we show why the variables g a b and f ab = 
\J—gg ab are preferable over g ab in choosing the dynamical variables for formulating an action principle for 
general relativity. Section 4 provides some useful relations between known expressions in relativity and 
our canonical variables. Section 5 discusses the variation of the surface term in the action and explores 
its relation with the canonical variables. Finally, we present our conclusions in Section 6. 

The conventions used in this paper are as follows. We use the metric signature (—,+,+,..,+). The 
fundamental constants G, h and c have been set to unity. The Latin indices run over all space-time 
indices while Greek indices will be used for purely spatial indices. The tensor density corresponding to 
a tensor will be denoted by the corresponding letter in calligraphic font. For example, the Lagrangian 
density will be denoted by C = y/^^gL. 




2 Preliminaries: The Structure of the Einstein-Hilbert Action 

In this section, we shall rapidly review some ideas and relations known in the literature, in order to set 
the stage. As we said, there is an inherent difficulty in trying to build an action principle for general 
relativity from which the gravitational equations of motion can be derived, which sets it apart from 
all other theories. The dynamical equations of motion are not expected to contain derivatives of the 
variable of order greater than two. This normally requires the action to contain not more than the first 
derivatives of the dynamical variables. But in general relativity, we would like the Lagrangian to be a 
scalar and any such scalar built out of the first derivatives will be zero in the local inertial frame, and 
by virtue of being a scalar will be zero in any other frame. The traditional way of dealing with this 
situation is to construct an action consisting of the metric and its first and second derivatives and then 
arrange matters such that the equations of motion are only second order in the derivatives of the metric. 
In fact, this demand alone is enough to uniquely identify the action in D = 4 dimensions. This action is 
the Einstein-Hilbert action, given by 



with L eh = R in our notation. It is useful to define a quantity Q a bcd and write the Lagrangian density 
in an equivalent form as: 





(1) 




Q, 



a 



bed 




(2) 
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The advantage of this form is that it is readily generalized to gravitational theories in more than 4 
dimensions where our conditions allow other terms in the action in addition to the Einstein-Hilbert term 
(see e.g., Chapter 15 in [32], [42]). It is a well-known result that £eh can be decomposed into a bulk 
term, which is quadratic in the derivatives of the metric, and a surface term which contains all the second 
derivatives (see [36,37]). The variation of the bulk term alone can furnish the equations of motion (as 
explicitly shown in e.g., chapter 6 of [32]) while the surface term, when integrated over a horizon, is 
related to the entropy of the horizon. The decomposition into bulk and surface terms is given by 

^EH = -Cquad + £sur, (3) 

where we have defined the bulk Lagrangian density and the surface Lagrangian density respectively 

as 

£ quad = 2j—gQ a bcd T a dk T k bc - C sm = 2d c [V=dQ a hcd T a bd ] . (4) 

We will also require the expression for the bulk Lagrangian in terms of the derivatives of the metric. 
This is given by: 

iquad = \M ab ^ k d a g bc d t g jk , (5) 
where we have defined a 6— indexed tensorial object M abci: > k , given by: 



M abcijk = gO-igbOgik _ \( g ai g bj g ck + g ai g cj g bk^ 



2 



2 

- \{g ak g bc g l} + g a3 g bc g lk + </'VV" + <?W 6 )- (6) 

This has the following symmetric properties. It is symmetric in b, c and j, k and also under exchange of 
the index triplets (abc) and (ijk). Eq. (5) allows us to obtain the canonical momenta corresponding to 

£*quad 

o(o a gbc) 2 

Note that M abc is not a tensor. Nevertheless, we can define a rule for raising and lowering of indices as 
in the case of tensors and lower the last two indices to define: 



Thus, \J—gM bc is the negative of the canonical momentum corespondent to g bc . The negative sign in 



L bc 

the last term arises because S(d a g bc ) = —g bd g ce 5(d a gde)- 

As we have stated already, the equations of motion can be derived from the bulk term alone. Once 
the equations of motion are obtained, we can find solutions to them, including black hole solutions, say, 
without ever bringing into discussion the surface term. If we now evaluate the surface term of the action, 
on the horizon, it will reproduce the entropy of the horizon in the Euclidean sector. (In general it gives 
tTS where r is the range of integration; we get S when r = /3, the inverse temperature.) The fact that the 
surface term, which is not supposed to know anything about the equations of motion, gives the entropy 
when integrated over a horizon demands an explanation. We stress that, this result — in fact — is a 
direct hint that gravitational action principle contains information about horizon thermodynamics. 

The algebraic answer to this question lies in realizing that the bulk term £ qua d and the surface term £ sur 
are not independent, but related to each other by the relation 



c g a b- 



d{d c g ab ) 



(9) 



Since it relates a quantity on the surface with a quantity in the bulk, this relation has been termed 
"holography" in the past [22-24]. It is the same —d(pg) structure of the surface term that allows us 
to derive the equations of motion by fixing the canonical momentum at the boundary rather than the 
metric [23]. We shall discuss this issue in the next section. 
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Using the canonical momenta M abc as defined in Eq. (7), we can rewrite Eq. (9) as 

C sm = - [d c (^gga b M cab )] = d c (V=gV c ) , (10) 

where we have defined a one-indexed non-tensorial object V c given by 

V c = -g ab M cab = g ik T? k - g ck V k n m = 2Q a bcd T a bd = 2Q^ kc d l9]k = ~- g d b {gg bc ). (11) 

Using the V c thus defined, we can write down the following expression for the canonical momentum 
M abc : 

M abc = g bd g ce T a de - \g bd g ac T de - \g cd g ab T de - \g bc V a . (12) 

Having thus set up the initial framework, we shall now turn to the discussion of how certain "holograph- 
ically conjugate" variables (HCVs) are more suited for the variation of the gravitational action. 



3 Holographically Conjugate Variables 

3.1 Variation of the Action and the Holographically Conjugate Variables 

Let us now consider the variation of the Einstein-Hilbcrt action. The variation of the Lagrangian in 
Eq. (3) is (see Chapter 6 in [32]) given by 

S£ EH = S({C quad } + {£ sur }) = {^—gG ab 5g ab - d c [^—gM c ab 5g ab ] } - 5{d c (^g ab M c ab )} 

= ^G ab Sg ab ~2d c [V^M c ab Sg ab }~d c [g ab S(V^M c ab )} (13) 

where M^ b is as defined in Eq. (8). Note that the surface variation has two parts, one part arising out of 
the variation of the metric while the other part contains the variation of y/—gM° b , which is the negative 
of the canonical momentum corresponding to g ab . Hence, in order to obtain the equations of motion 
we have to fix both g ab and its corresponding canonical momentum at the boundaries. But as we have 
argued in the introduction, there are some fundamental difficulties with this program. 

On the other hand, if we write the variation in terms of g ab , then a crucial sign flip leads to the cancellation 
of two terms (which added together in the previous case) and we find 

S(y/=gR) = -V=gG ab 5g ab - d c [g ik 6{^M cik )], (14) 

Thus we only need to fix the canonical momenta ^—gM clk , corresponding to gi k , at the boundary to 
obtain the equations of motion. So, some of the concerns raised in the introduction about the variational 
principle in general relativity can be addressed by treating the action principle as a momentum space 
action and using the covariant components of the metric, g ab as dynamical variables. [The addition of 
the term —d(pq) has a direct interpretation in quantum theory [22]. Usually, the propagator G{q2 1 q\) 
for the dynamical variable q can be obtained from a path integral using the action built from a quadratic 
Lagrangian L q (q, dq). The propagator G{p2,Pi) in momentum representation can be similarly obtained 
from a path integral using the action built from L q {q 1 dq) — d(pq).] 

Since we know that M clk is made of the metric and its first derivatives, it is natural to ask how the 
surface variation term in Eq. (14) looks like in terms of variations of the metric and the affine connection, 
in the spirit of Palatini. (We shall discuss the Palatini variational principle in Section 4.4.) The required 
expression is 

- d c [g, k S(V^gM" k )} = -d c [2^—gg bk 5{Qt d e Y dk )] = -d c [^—gg bk 5N^ k ] , (15) 

where Q ab = ^(8 a 5 b — S d S b ) is obtained by lowering an index in Q l a cd as defined in Eq. (2). We have 
introduced here a new object, 

n& = -n c + \ (itA Q + i1a q ) - Q a A + Q a Jn d , (is) 

which is made purely from the affine connection and does not involve the metric. The reason for 
introducing this object is not just aesthetics as will be clear later, in Section 3.3, along with the reason 
for the demand for symmetrisation in the lower indices. 
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Let us compare Eq. (15) with Eq. (14) and Eq. (13). We obtained some level of simplification in going 
from Eq. (13) to Eq. (14). But y^—gM cab , in spite of its simple interpretation as the canonical momen- 
tum, is a complicated object in terms of the metric and its derivatives (Eq. (7)) or in terms of the metric 
and the afhne connection (Eq. (12)). Hence, it is a pleasant surprise to find that the surface variation 
turns out to be a pure variation of the affine connection. As a consequence, we can let the metric be 
arbitrary at the boundary and fix just the connection at the boundary to obtain the equations of motion. 
If we consider the metric and the afhne connection as two different fields a la Palatini, then this is a 
perfectly acceptable situation, although a little peculiar since we need to fix only one field out of the 
two that we are considering. The reason for this can be traced to the fact that the Einstcin-Hilbcrt 
Lagrangian can be expressed in terms of the metric, the affine connection and the derivatives of the 
affine connection without involving the derivatives of the metric. 

If one wants to think of the variation in terms of the metric and its derivatives instead, defining the 
affine connection a priori in terms of the metric and its derivatives, as we want to do, then it is more 
difficult to justify setting the surface variation to zero at the boundary. As regards the possibility of 
imposed boundary conditions being incompatible with the equations of motion, note that the momenta 
have the same number of components as the derivatives of the metric and hence we should be able to 
leave the metric components as arbitrary and achieve the momenta constraints just by manipulating the 
derivatives. Thus, that issue seems to have been addressed, at least at first sight. The second issue, 
related to the uncertainty principle, is not applicable here as we have to fix just the momenta among the 
canonically conjugate variables. 

Thus, we can make a case for the use of g a b as the dynamical variables as opposed to g ab . We shall call 
the pair (s (I ^,/ = JM m ' , ) as "holographically conjugate" variables (HCVs) in the spirit of the —d(pq) 
structure of the surface term in terms of these variables. In the next section, we will see that the 
same conclusion can be arrived in an simpler manner using some scaling arguments. In the subsection 
after that, these arguments will lead us to another pair of "holographically conjugate" variables with a 
symmetric, contravariant, two-indexed object taking the place of g a b- 



3.2 Holographically Conjugate Variables through Scaling Arguments 

In this section, we shall show that the results of the previous section can also be obtained using simple 
scaling arguments. There is an alternate method of arriving at Eq. (9) using scaling arguments (see 
Project 6.3 in [32]). We shall describe this method below since it elucidates the role played by the use 
of g a b as the dynamical variable rather than g ab . 

Consider a Lagrangian L(qA,diqA), where qa is some dynamical variable and A denotes a collection of 
indices, such that the Lagrangian is a homogeneous function of degree \i in qA and degree A in diqA- The 
Euler-Lagrange function (see Section 4.3) obtained from L for the variable qA is 



P A _ 



dL 



0; 



dL 



d{diq A ) 



(17) 



Forming the contraction qAF A and using Euler's theorem for homogeneous functions, one can easily 
show that 



q A F A = (A + y)L - d t 



dL 



<1A 



d{diq A ) 



If we take L to be y 



-gLqund and qA to be g a b, we would have = — 1 and A = 
= —y/~ r g~G ab and qAF A = y^gR. Hence, Eq. (17) becomes 



rr gR = V^ffiquad - d c g ab 



gLqu&d \ 



d(d c g a b) J 



On the other hand, for the case of g ab as the variable, we have /i 
qAF A = — ^^gR. Thus, we obtain 



-3, A 



2. F 1 



-gR 



"5-kquad + 8 C ( g 



d(d c g ab ) J 



-ff-kquad + 9 C (g a bV-9 M 



cab\ 



(18) 

Also, F A = 

(19) 
(20) 

v ^gG ab and 
(21) 



G 



with a crucial sign difference in the second term. 



In both cases, the second derivatives are confined to the surface term. Hence, if the dynamical variables 
and their derivatives are fixed on the boundary of the volume under consideration, the equations of 
motion can be obtained by varying £ qua d alone and will be of second order in the derivatives of the 
metric. As has been explained in the introduction, we cannot fix both the dynamical variable and its 
corresponding momentum at the boundaries. Now, the surface term that arises from the variation of 
£ qua d alone will be of the form d(pSq) and can be eliminated by fixing the dynamical variable q at the 
boundary. But £ qU ad has the disadvantage that it is not a tensor density and, in fact, vanishes in a local 
inertial frame. A —d(pq) term added to this Lagrangian density, as in the case of q = g ab in Eq. (19), 
will make it into a tensor density and also modify the surface term in the variation to the form —d(q5p), 
allowing us to fix just the canonical momenta at the boundary and obtain the equations of motion (see 
section II in [23] for a detailed discussion). The explicit variation of Eq. (20) is 

S(^gR) = -^~gG ab 8g ab - d c {g lk 6(V=gM ctk )}, (22) 

showing that we only need to fix y/—gAf clk at the boundaries to obtain the equations of motion. This 
is equivalent to fixing the variation of the connection at the boundaries, as can be seen from Eq. (15). 
This insight suggests that Hilbcrt action is better considered as describing a theory in the space of the 
canonical momenta corresponding to g ab . 

On the other hand, Eq. (21) tells us that the surface term occurs with the wrong sign when the variable 
g ab is used. The explicit variation in this case leads to 

5{J—gR) = V^G ab Sg ab - d c \g ik 6(V=9~Mf k )) - 2d c [^Mf m 8g lm ). (23) 

The surface term now has the variations of both the dynamical variables and the momenta. Generally, 
textbooks choose g ab as the dynamical variable because some extra care has to be taken to derive the 
equations of motion if g ab is taken as the variable (see Exercise 6.9 in [32]). But then one would either 
have to ignore the variation in the surface term or cancel it with a counter-term like the Gibbons- 
Hawking- York counter-term [31,34], but neither option is as simple and neat as using g ab as the variable 
and fixing the canonical momenta at the boundaries. 

In the next section, we present another object which, along with its canonical momentum, forms a pair 
of HCVs (holographically conjugate variables) and can be used in place of g ab for a variational approach 
to general relativity. 



3.3 An Alternate Pair of Holographically Conjugate Variables 

A natural question that arises is whether g ab is unique in providing us with a neat and clean variational 
principle for general relativity. In Appendix A, we have used scaling arguments to investigate if there 
are other variables which may be used in place of g ab in the variational approach to general relativity. 
Although g ab cannot be considered as a good variable for this purpose, we have found an object with 
two contravariant indices which appears suitable. This object is the tensor density \J—gg al ', denoted in 
this paper by f ab for typographical convenience. The variable f ab scales linearly with g ab i.e g ab — > ag ab 
leads to f ab —> af ab . It is precisely this linear scaling that makes f ab a suitable substitute for g ab , as 
demonstrated in Appendix A. We shall use the symbol f ab for the corresponding covariant tensor density, 
V~99ab- Note that f ab is not the inverse of f ab . See Appendix B for some useful properties and relations 
pertaining to f ab and its variation. 

As proved in Appendix A, writing the Einstcin-Hilbert action in terms of the new variable f ab reproduces 
the structure of Eq. (19). We obtain 

V~gR = V^iquad - dcir ^^jlf ] = V^fliquad - d c [f ab K b ] , (24) 

where we have defined the object Nf k to be the canonical momentum corresponding to f ab for £ qua d 
by: 

W - ^(y^ffZ/quad) , s 

Nik = d(d c f») ■ (25) 
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The explicit expression for Nf k has already been written down in Eq. (16). The reason for the demand of 
symmetrisation in Eq. (16) is evident now. Note that it is a simpler object than M clk and is constructed 
from the affine connection alone. 

Next, let us look at the expressions for variation of the Einstcin-Hilbcrt action in terms of f ab and N£ b . 
As we have argued in Section 3.2, the —d(qp) structure in Eq. (24) gives rise to a ~d(q5p) surface term, 
along with the bulk term which provides the equations of motion as usual. Using Eq. (14), Eq. (15) and 
Eq. (B.6), we can write down the variation as 

5(V=gR) = Ra b Sf ab + f ab SR ab = R ab 5r b -d c [r k 5N^}, (26) 

In addition to their utility in simplifying the variation of the action, there are two main advantages to the 
HCVs (holographically conjugate variable pairs) (g a b, \/~gM cab ) and (f ab ,N£ b ). The first one is that 
many known expressions and formulae simplify considerably when written in terms of these variables and 
make our theoretical life easier. The pair of variables (f ab ,N^ b ) gives a particularly stellar performance 
in this regard and will be our variables of choice for most of the work in this paper. More importantly, the 
variations of these variables on a horizon will be shown to have a direct thermodynamical interpretation. 
In the next section, we shall demonstrate the first point by writing down several well-known expressions 
and formulae in terms of the HCVs (holographically conjugate variables pairs). 

Historically, it was indeed noted in the early days of general relativity that f ab is a good variable to use 
and was even given preference over g ab at times [36,37]. The variables f ab , N£ b were later used, albeit the 
non-symmetrised version of N£ b , by Einstein in his work with Kauffman attempting to go beyond general 
relativity [39,43]. A recent paper that uses these variables but does not make the identification of N£ b as 
the canonical momentum corresponding to f ab is [44] . The paper [45] contains the analogues of some of 
the expressions in terms of these variables given in the next section, but the calculations have been done 
as perturbations on a flat metric. We believe the holographic relation between the bulk and the surface 
terms of the action gives a simple and elegant motivation for using the variables (f ab ,N^ b ). 



4 Gravity in terms of f ab and N£ b 



In this section, we shall show that many objects of common use have simpler expressions in terms of f ab 
and N^ b than in the conventional description. First, note that Eq. (B.4) allows us to write 

dj ab = V=gB^d c g lm ; B% = \{5 l a 8T + 8 l b 5™) - (l/2)g lm g ab (27) 

and since Bf b n is independent of the derivatives of g ab , we can write 

d{dcfab) =tiV=gBt>. (28) 



d(d d g lm ) 

Using this relation and Eq. (8) and Eq. (25), we get the relation 

V=gM c ab = -V=gB%N? m , (29) 

which enables us to relate the canonical momenta corresponding to g ab and f ab by: 

M cab = _ B ab N cde ( 3 q) 

This can be easily inverted using Eq. (B.3) to give 



jy-cab 



= -Bf e M cde = -{M cab + ^g ab V c ). (31) 



Looking at Eq. (12), we can see that N bc is related to the Christoffel symbols by the simple expression 
Eq. (16). Thus, we can infer that although N bc is not a tensor, 5N£ b is a tensor because ST^ is a tensor. 
Eq. (31) can be readily converted to an expression in terms of the first derivatives of the metric using 
Eq. (7). This expression is 

N cab = —-Bf e M cdeijk dig jk . (32) 



8 



We can also write N£ b in terms of d c f ab using the inverse of Eq. (27), to obtain: 

N c ab = -=L [g ci B abrs - 5 c r (5 l a g bs + S l b 9as)] dj rs (33) 

To replace the first derivatives in the theory with canonical momenta, we need to invert Eq. (16). This 
is easily done by assuming the following ansatz for the Christoffcl symbols: 

Kb = *K b + b{N d ad 5 c b + N* d 5 c a ). (34) 
and substituting back in Eq. (16) to solve for a and b. We obtain 

Kb = -Kb + \{N d ad 5l + Nfal). (35) 

Now we can substitute for N£ b in terms of M^. (We shall not display this equation explicitly as it 
does not appear to simplify further.) We can use the relation d c g a b = Ta.&c + {a f-> 6} and obtain the 
derivatives of the metric to be: 

d c gab = {-g ae {Nt c - \{Nt d 5 e b + Ng&)]} + {a^ b}, (36) 

and, further 

gjab = [f ad (N b c _ l _ 5 b N C de)] + [a ^ b] _ (37) 

We shall now make use of these expressions to express the Lagrangian, curvature tensor etc. in terms of 
the HCVs (holographically conjugate variables). 



Ra 



b 



4.1 Riemann Tensor, Ricci Tensor and Ricci Scalar 

We next give the formulas for the Riemann tensor, the Ricci tensor and the Ricci scalar in terms of N bc 
for ready reference: 

RU = [-(Std c + N b e J(NS e -^S a d Nl f ) + ^S a c N§ e N^ f ]-[c^d\ (38) 

- (d c N c ab + N a ad Ni - \K c Ni d ) (39) 
R = -g ab d c NC b -L quM (40) 

v^i? = -r h d c K b - l -K b d c r b (4i) 

= \Kbdcf ab -d c (f ab Kb) (42) 
= ~[f ab doK b + d c (r b N c ab )] (43) 

Here, Eq. (42) is the usual decomposition of ^—gR into a bulk term and a surface term as given in 
Eq. (3) while Eq. (43) is an alternate decomposition in which the bulk term contains second derivatives 
of the metric. 



4.2 The Lagrangian 

We can substitute for Ts in the bulk term in Eq. (4) in terms of N£ b s from Eq. (35) and obtain the bulk 
part of the Lagrangian to be: 

i quad = g bd N dj N* [8\8{ - i^] = g bd [rr[N b N d ] l -Tr[N b ]Tr[N d ]^ (44) 

where we have taken N£ b to be the cb th element of a matrix N a , and Tr[N a ] is the trace of this matrix. 
This also allows us to write: 

£quad = \K b d c f ab , (45) 
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with striking simplicity. Note that this exhibits a 'p<?/2' structure, which is a consequence of the fact 
that the Lagrangian is quadratic in q. We can also write L qua< i in terms of M£ b as 



V.ad = 9 hd M c dk Mi - 2 -g ir MtM? lk - ±g bk M l bl M d k + l -g dk g ir M? r g*« M* y 

= g bd Tr[M b M d ] + ^Tr[A/ fc ]l/ fc - \g bd Tr[M b ]Tr[M d ] + \g bd V b V d , 
3 3 6 



(46) 



which, in comparison to Eq. (44), appears to be quite complicated. But, from Eq. (5) and Eq. (7), the 
l pq/2' structure arises with the variables (g a b, \/—gM cab ) also, as given below: 



quad 



-9 



M cab d c g ab 



The formulas for L sur have already been given in Eq. (10) and Eq. (24). We have 

£sur = d c {J—gV c ) = -d c [g ab (V^M cab )] = ~d c [.f ab N c ab )]. 
In fact, the above relations hold even without the derivatives, 

V~gV c = -g ab {J—gM cab ) = -f ab K b , 
This equation will be of use in the study of horizon thermodynamics. 



(47) 

(48) 
(49) 



4.3 Euler Derivative and the Equations of Motion 

The Euler derivative of any function K\<p, di4>, ...] with respect to the variable <f> is defined as 



SKfadit,...] _ dKfadit,.. 



O0 



d a 



d(d a( j>) 



d a db 



d{d a d b 4>) 



(50) 



The function thus obtained is also called the Euler-Lagrange function resulting from K for the variable 
(f>. We will use the notation E[K, <f>] for this object. We have already made use of the Euler-Lagrange 
functions in Section 3.2. 



The most general variation of K can be written as 

surface variations 



SK = ^ 



(51) 



Thus, when the surface terms can be consistently put to zero, the equations of motion are obtained 
by equating the Euler derivative to zero (which is the origin of the nomenclature). We shall now list 
the Euler-Lagrange functions obtained from £ qua d for the dynamical variables under our consideration, 
namely g ab , g ab and f ab . The functions are, respectively, 



E[£ qu3 ,d, g ab ] — \f—gG ab , 
E[Cq Uad ,g a b] = -V^gG ab and 

-E[£quad, f ab ] = Rab ■ 



(52) 
(53) 
(54) 



The equations of motion in the absence of matter can be obtained by equating these functions to zero. 
In the presence of matter, the matter Lagrangian has to be added to £ qua d and the Euler-Lagrange 
function of the total Lagrangian has to be obtained. We then get the equivalent expressions: 



Gab = 87rT Qb ; G ab = 8nT a b] R a b = 87r(T ab — -^T) 
From Eq. (39), we can write down the equations of motion as 

d c N c ah = -N^Ni + ^NlMd-MTab-^T) , 



(55) 



(56) 



a first order differential equation in N£ b with a matter source term. In the local incrtial frame, the 
equation becomes 

d c N ab = -8n(T ab - Mt) (57) 
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which gives a conservation law for the canonical momenta N£ b valid in a space-time volume small enough 
to respect our local inertial frame. The corresponding equations of motion in terms of M cab look even 
nicer in the local inertial frame. They are G ab = d c M cab = %nT ab i.e. 

d c M cab = 8irT ab (58) 

Note that in the local inertial frame y/—g is unity and hence M cab = ^/—gM cab plays the role of 
canonical momentum corresponding to the metric. Thus, if we take a region small enough for a local 
inertial frame to be enforced, the surface integral of the canonical momenta corresponding to the metric 
components gives the components of the matter energy-momentum tensor contained in that volume. 
Such a balance between gravitational variables and matter variables will again arise when we study 
horizon thermodynamics. 



4.4 Palatini Variational Principle 

In classical mechanics, the variation of the action is carried out regarding the dynamical variable q and 
its first time-derivative, q as the independent variables, for the class of actions which do not depend 
on the higher time-derivatives of q. The momentum p is then defined as the partial derivative of the 
Lagrangian with respect to q. There is an alternate way of carrying out the variation considering q and 
p as independent variables, with the Lagrangian being defined as L — pq — H(q,p), where H(q,p) is the 
Hamiltonian of the system [46,47]. In this case, after fixing the variations of q at the boundary, the 
variation with respect to p gives the equation q = dH/dp, while the variation with respect to q gives 
p = —dH/dq, which are the well-known Hamilton's equations. This variational principle is referred to 
as modified Hamilton's principle. 

Analogously there are two well-known variational principles in general relativity: one in which the 
components of the affine connection are considered as given in terms of the derivatives of the metric and 
the variation of the Einstcin-Hilbcrt action is carried out in terms of the variation of the metric and its 
first and second derivatives; and one in which the metric and the affine connection arc considered as 
independent and varied separately. The second method is called the Palatini variational principle [48]. 
But, as we have been arguing the case for the use of the HCVs (holographically conjugate variables), 
we shall now outline a version of the Palatini variational principle in terms of the variables (f ab ,N£ b ). 
As these form a (q,p) pair, this would be a direct analogue of our classical mechanics exposition of the 
alternate variational principle. This aspect is in contrast with the usual approach in which the metric 
and connection are treated as independent variables, because the connection T bc is not the canonically 
conjugate variable to the metric g bc . 

Using Eq. (39), the Einstein-Hilbert Lagrangian can be expressed as 

J—gR = f ab R ab = f ab {-d c N c ab - KdNi + ^N c ac N^ d ) . (59) 

We will now vary this Lagrangian considering f ab and N l ab as independent variables. Note that R ab is 
a function only of N l ab and is independent of f ab , which should not be surprising as we know R ab to be 
a function only of the affine connection. The variation of the action (Eq. (59)) with respect to N l ab is 
given by 

cab 



gR)\ f *> = f ab 5R ab (60) 
dcfab _ 2f ad N b d + lf*™ N $ mS ^ 5K b - d c {f ab 5Nl b ) . (61) 

Once we fix N bc at the boundary, we obtain the corresponding equations of motion by equating the 
symmetrised coefficient of 5N£ b to zero. These equations are 

gjab = f ad Nbcd + fb d N a d _^ f a mN d m6 b_^ fbmN d m§ a ( g 2) 

Note that Eq. (62) is identical with Eq. (37). Thus, the action principle dictates the connection between 
N£ b and (f ab , d c f ab ). In order to connect this up with the standard result V ' c g ab = obtained from Pala- 
tini variational principle in its conventional form, we need to keep in mind that the covariant derivative 
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is, as of now, defined as usual in terms of the affine connection but the afhne connection is not related 
to the metric or its derivatives. Substituting f ab — y/—gg ab in Eq. (62) and contracting both sides with 
g a b, we can obtain the relation N d d = {3/4)g ab d c g a b which gives us: 

V~gN* d = ^d c V~g; i.e., V~g^ d cd = d c V~g (63) 

We shall now use this result to evaluate yj— gV c g ab . Expanding yj— g'\7 c g ab , we find 

^V,g ab = d c f ab - g ah d c y[—g + f mb T a cm + f ma T b cm . (64) 

Next, using the identity dcyf 1 ~g = (2/3)y/^gN d d from Eq. (63) and expressing the connections in terms 
of N ab variables by Eq. (16), we obtain 

sf—gV c g ab = d c f ah f ad N b d - f bd N a cd + l -f am N d md 8 b + ^f bm N d nd S a c , (65) 

which vanishes due to Eq. (62). Hence, the equation of motion (Eq. (62)) is precisely the metricity 
condition obtained by the variation of the connection in conventional Palatini formalism. 

We shall now provide an alternate perspective on these results which would prove quite fruitful. Note 
that Eq. (62) can be written as 

d cf ab = - JLr (V=gR + f ab dcK b ) (66) 



9N c ab 

We have obtained here an analogue of the Hamilton's equation q = dH /dp. The analogy can be made 
more precise as follows. Let us define: 

n g = f ab (K d Ni - lK c N d d ) (67) 

Then, with the notional correspondence f ab — > q and N£ b — > p, we can establish: 

yf^gR — > -qdp -U g = {pdq - U g ] - d(qp) = £ quad + £ sur (68) 
Comparing Eq. (67) with Eq. (44), we see that 

Ug = £quad > T^pdq (69) 

Thus, the quadratic Lagrangian density that is used to derive the equations of motion can be also be 
interpreted as a Hamiltonian density. The equation Eq. (66) can then be rewritten in the desired form 
as 

d c fo,b = (70) 
dN c b V ) 

Proceeding by analogy, our next stop would be the other Hamilton's equation of motion, p = —dH/dq. 
Consider the variation of the action (Eq. (59)) with respect to the q variable f ab , given by 

5(y=gR)\ Nib =R ab 6f ab . (71) 

The equation of motion obtained on extremising the action with respect to variations in f ab is R ab = 
which is the same as: 

- K d Ni + \N c ac N d d = d c N c ab (72) 
Referring back to Eq. (67) , we see that this equation can be expressed as 

9oK b = (73) 

giving the second of the Hamilton's equations. But unlike the case in classical mechanics where the 
momentum would be conserved in the absence of external forces, we see that N^ b is capable of driving 
its own change, due to the nonlinear nature of gravity. 
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The next natural step is to consider the inclusion of the matter Lagrangian density £ m = yf—gL^. This 
can be accomplished by defining a total Hamiltonian as 



Ra — Cr, 



(74) 



If the first term R g , which is equal to £ qU adj be considered as a kinetic term, then it is natural to think 
of £ m as a potential term as far as gravity is concerned. We shall make here the assumption that the 
C m under consideration depends only on f ab and not on N% b . In such a case, Eq. (62) retains its form. 
Thus, if we choose to express everything in terms of the metric and its derivative, our assumption is 
tantamount to the assumption that C m does not depend on the derivatives of the metric. (This is similar 
to the case in classical mechanics when we consider velocity- independent potentials.) 



If we now we take the usual definition of the matter energy-momentum tensor as 

2 8(^L m ) 



g 8g ab 



we can obtain the following equality 

<9£ m 



Of 



ab 



; T 9ab 



jDlrnrp 



(75) 



(76) 



where B a ™ was defined in Appendix B. We have defined a new object T ab here, which bears to T ab the 
same relation as G a b bears to R ab . 



Eq. (73) is now modified to read 

ft M c - 
°c ly ab- gj-ab 



d(167r£ n 



df 



ab 



= -K d Nt + -N c ac N« d - 8irT ab 



(77) 



Thus, the presence of matter introduces an extra source term in the equations governing the evolution 
of N£ b - Using Eq. (39), it is easy to verify that Eq. (77) is equivalent to the usual Einstein's equation of 



motion, R ab = 8nG[T ab - (l/2)Tg a 
Finally, note that the variation of the total Lagrangian, gravitational plus matter, on varying f ab and 

N ab is s iven b y 



6{y/=gR+ 167r£ r 



gjab _ fadN b d + l f 

-d c (f ab SN c ab ) 



am \jd cfc 
N dm°c 



5N c ab - {(d c N c ab + N c ad Ni - \N c ac Nt d ) + 87rT afc ] ( 



(78) 



The surface variation given in the second line of the above equation contains only variations of N£ b and 
not of f ab , which happened essentially because the Lagrangian with which we started with does not 
have derivatives of f ab . Thus, we need to fix only the "momenta" N£ b at the boundary to obtain the 
equations of motion. 



4.5 Noether Current 

So far, we have introduced two sets of naturally conjugate variables and expressed some relevant quanti- 
ties in general relativity in terms of them. The usual Einstein's equations have also been derived from the 
variation of these variables. In this subsection, we shall relate the Noether current — which arises from 
the diffcomorphism invariance of the action — with the variations of N bc and f ab . For a general covariant 
Lagrangian L(g ab , R a bcd), the Noether current, corresponding to the diffcomorphism x a — » x a + £ a , can 
be shown to be related to the Lie derivative of F^ c with a particular contraction of the indices with the 
quantity P abcd = dL/dR abc d- Since N bc is a linear combination of the connections, we can easily convert 
this relation into the required relation. 

To do this, we first recall the general form of the Noether current for a Lagrangian L(g abl R a bcd), corre- 
sponding to the diffcomorphism x a -> x a + £ a . It is given by (see, e.g., Project 8.1 of [32]). 

J a = 2E^ b + LC + 5 iV a , (79) 
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where 8^v a is such that V a {8^v a ) is the surface term in the variation of the Lagrangian under the 
diffcomorphism and E b — is the equation of motion. For the Lanczos-Lovelock models of gravity, for 
example, E e b l = P aki ^R bkij - \8%L and 8 ( v j = 2P lb ^ d \' b S i g dl (see Section 15.4 in [32]), where P abcd = 
(dL/dRabcd) as already mentioned. The tensor p abcd has the algebraic properties of the curvature tensor 
and, additionally, it is assumed to be divergence-free in Lanczos-Lovelock models: \7 a P abcd = 0. Hence, 
for Lanczos-Lovelock models, the relation between the surface contribution of the variation of action and 
the Noether current is 

S^v a = J a - 2Tl a h i h , (80) 

where 11% = /'" , *'< ///,;,, ,. We now switch to the language of Lie derivatives. For an indexed object A, 
while 8% A = A'{x) — A(x) is the functional change in A at the same value of the space-time coordinates 
under an infinitesimal coordinate transformation from x % to x l + , the Lie derivative of A is defined as 
£^A = A(x) — A'(x), thus giving = —8%. Now, we can use £^g d i = Vd£, + V,^d to find 

£^v a = -2P abdl V b V d ^ - AP aibd V b V d ^ = -2P abdi V b V d & + 2R a m C n , (81) 

where, in the last step, we have used P ai6rf V 6 V d & = -(l/2)P albd R mlbd £, m . Therefore, the Noether 
current turns out to be 

J a = 2Kl£ b + S iV a = 2TZ^ b - £ t v a = 2P abdi W b W d £, l , (82) 

We will now re-express this in terms of the variation of the connection. We know that, although T bc is 
not a tensor, its Lie derivative is a tensor: 

£^1 = V b V c C + R a cmb C n . (83) 
Contracting the above expression with P lbc a we find 

2P^ a £^ a bc = J 1 - 27Cr ■ (84) 
Therefore, from Eq. (80) and Eq. (84), 

Sp* = 2P lbc a £^ a bc . (85) 



It turns out that, in general relativity, we can easily write this expression in terms of N bc . To do this, 
note that in the case of general relativity, L = R and P lbca = \{g lc g ab - g la g bc ) = Q lbca . Now, since N£ c 
is related to T bc by Eq. (16), it is easy to show that 



2P ibc a £^t c = g xy £ s Ni y . (86) 

So, the surface term in the variation of the Lagrangian density y/~gR can be expressed in different ways 
as (compare with Eq. (15)) 



-gV^v 1 ) = 2^V i (P^ a £ i T a bc ) = x — V,. :.-,'•''/', .Y; v ) , (87) 
while the sought-after expression for the Noether current in terms of the variation of N bc is given by 

J' =21Z l m C l + g xy £ i K y (88) 



4.6 Canonical Energy-Momentum Pseudotensor 



From the bulk part of the Lagrangian, one can define the canonical energy-momentum pseudotensor (up 
to overall factors, also known as the Einstein pseudotensor [36,37,49,50]) as: 



d{y/—gL quad ) ■ 
— — \ — cikgab - <V 



d(dj ab ) 



-d k f ab 



f -Lquad 



N ab d k f ab - -si(K b dj ab ) 7 



(89) 
(90) 
(91) 
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which, when taken together with the matter energy-momentum tensor, satisfies the conservation law (see 
Appendix C for a proof): 

$(4 - 16ttV=M = 0. (92) 

Finally, we mention that the "trace" of the pseudotensor can be related to the Hamiltonian T-L g , defined 
by Eq. (67). From Eq. (91), after contracting with 5 k , we find t\ = —N^ b d c f ab which is — 2£ qua( j (see, 
Eq. (45)). Therefore, using Eq. (69), 

4 = -2n g . (93) 

This concludes our discussion of standard features of general theory of relativity in terms of the canoni- 
cally conjugate variables (f ab , N£ b ). We shall now discuss the connection between the HCVs (holograph- 
ically conjugate variables) and the thermodynamics of horizons. 



5 Thermodynamics with the Holographically Conjugate Vari- 
ables 

In this section, we shall show that the variations p5q and qSp obtained from our sets of HCVs (holograph- 
ically conjugate variables) have direct thermodynamical interpretations when integrated over horizons. 
We shall first prove the results in a general static spacetime. Then, we shall specialize to the spherically 
symmetric case and examine Schwarzchild and Rcissncr- Nordstrom horizons in order to obtain a physical 
feel of our results. Finally, we shall show how our results generalize to integrals over any null surface 
which acts as a local Rindlcr horizon. 



5.1 Preliminaries 

We shall first set up a couple of relations involving the two canonical momenta before we proceed to the 
thermodynamic relations. From Eq. (22), Eq. (26) and Eq. (B.6), we see that 

d c [g lk S(V^M^ k )] = d c [J—gg bk 5N c bk ]. (94) 

Since V c = -g ab (yf=g~M cab ) = -f abN ab> wc also have 

8g lk {V=gM clk ) = 5f ab K b . (95) 

The variations 5(^—gM cxy ) and g bx g ky 6N£ b are not equal but become equal on contraction with g xy . 
To characterize the difference, we define a tensorial quantity 

H ab,cd = g ac g bd _ Q/Q g ab g cd (96) 

which has the property that H ab ' cd g a b = H ab,cd g c d = 0. This quantity is also a projector since we have 
f{abffcd,ef _ jjab.ef j n an y indexed quantity Q ab - can be written as 

Q ab - = \g ah g cd Q cd - + H« b Q cd -. (97) 



where dots indicate indices which are not displayed. A contraction with g a b will pick up a contribution 
from 
write 



from the first term and a contraction with HH will catch the second term. Using H ab,cd , wc can 



5(y^M cxv ) - [g bx g ky 8N^ b ] = 2(H ik ' xy 6 c z - H lk > cy 8 x )6T z xy , 

which makes it clear that S(y/—gM cxy ) and g bx g ky 5N^ b are not equal but become equal on contraction 
with g xy . 

Having set up these relations which would allow us to easily hop between the two pairs of HCVs (holo- 
graphically conjugate variables) that wc have, we shall now study the variations of the surface term of 
the Einstcin-Hilbcrt action on a horizon in terms of our variables. 
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5.2 Surface Term and its Variation 



In this section, we shall calculate the variation of the surface term for an infinitesimal change of the 
metric. It will have two terms: one corresponding to the variation of the metric variable - like pSq; while 
the other one is like q5p, coming from the variation of the momentum variable N£ b or y/—gM cab . We 
first give the general expressions for these terms in terms of metric coefficients and the perturbation h ab . 
Next, we shall explicitly calculate them on the horizon for a general static spacetime. Interestingly, the 
first term will lead to TSS while the other will give S5T, where T and S are the temperature and entropy, 
respectively, of the horizon. Finally, we will generalise this result to an arbitrary null surface. 

As already described in Section 2, the surface term in the Einstein-Hilbert Lagrangian is given by C sur = 
d c (^J—gV c ), where V c is defined in Eq. (11). In order to look at the variations of this surface term, we 
need to compute 5(y/—gV c ). From Eq. (49) and Eq. (49), we can split the variation of the surface term 
into two components as 



-NZ b 5f ab - 5N c ab r» 
-^M cab 8g ab - 5(y/=g~M cab )g ab 



(99) 
(100) 



In fact, from Eq. (94) and Eq. (95), we know that N r ab 5f ab = y/=gM rab Sg ab and 5N r ab f ab = 5(y/=gM rab )g ab . 
We shall work with the variables / and N£ b for the moment. In terms of these variables, it is clear 
that the first term in the above set of equations corresponds to variations of the metric while the second 
term corresponds to variations of the connection. 

To facilitate the calculation, we shall write down expressions, up to linear order, for small changes in the 
metric. Suppose the change in the metric is of the form g ab — > g ab + h ab with h ab treated as a first order 
perturbation. (To be more precise, we should work with eh ab , and retain terms linear in e and finally 
set e = 1. We shall not bother to do this.) Under this change, the terms in Eq. (99) arc evaluated up to 
linear order in h ab as 



and 



SN%f k 



N 

P k 
1 



f k 6f k = V=9N? k (±hgi k 



(101) 



9 ab djh bk + \g ah d b h dk + h ab d ]9bk - hi ab d b9ok 



-J 



ah 



g mn d k K 



l d k9n 



(102) 



where h = g ab h ab . In the following, we shall evaluate the above terms on the horizon for a general static 
spacetime. 



5.2.1 A general static spacetime 

For any static spacetime with a horizon, we can choose an arbitrary 2— surface and write the line element 
in the form [51]: 

ds 2 = -N 2 dt 2 + dn 2 + (T AB dy A dy B , (103) 

where n represents the spatial direction normal to the specified 2-surfacc and o ab is the transverse 
metric on the 2-surface. We shall assume that there exists a Killing horizon determined by the timelikc 
Killing vector £ = dt with the location of the horizon given by the condition A^ 2 — > 0. We will choose 
the coordinates such that n = on the horizon. Then, near the horizon, N and aAB has the expansion 
[51], 

N = nn + G(n 3 ); a ab = \o H {v)\ab + ^My)Us" 2 + 0(n 3 ) . (104) 

To evaluate Eq. (101) and Eq. (102) on the horizon, we shall first calculate them on the n=constant 
surface and then take the limit of n — > 0. The non-zero components of h ab are h u = —2N5N; Hab = 
5a ab and the relevant non-zero components of N bc are 

NH = -Nd n N; Kn = ^+ l -a AB d n a AB ; N n AB = \d n a ab ■ (105) 
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Using all these in Eq. (101) and Eq. (102) we find 



N? k 5f jk = ^a AB 5a AB d n N -^^a Ac a BD 5a CD d n o AB 

+ Y^ AB Sa AB a CD d n a CD + V^a AB d n a AB SN ; (106) 



(107) 



aim 

P k 5N? k = 2^d n (SN)+N^a AB d n (8a AB ) - ^^o AC o BD 8o C Dd n o AB . 

Let us now evaluate Eq. (106) and Eq. (107), using Eq. (104), near the horizon. This will lead to 

N? k 5f* k \ H = 2k8{^); f k 8N? k \ H = 2^8K. (108) 

Integrating over the transverse variables and introducing the appropriate numerical factor, we find that 
the variation of the surface term in the Einstein- Hilbert action Aeh is given by (see Eq. (1), Eq. (3) and 



Eq. (24)) 



-L [d?x ± N: b 5r h = ±S(^)=TSS; (109) 

1D7T / 27T 4 



^JdWr b S Kb = (^)S(^) = SST. 



(110) 



In obtaining these results, we have appealed to the zeroth law of black hole thermodynamics [52,53] (Sec 
also Eq.(61) in [51]) and taken k and Sk to be independent of the transverse variables. 

Before proceeding further, we address, as an aside, the following question: We saw in Eq. (110) that 
there is a clear correspondence between the variations p8q and qSp on the one hand and TSS and SST 
on the other hand. How special are the variables (g a b, y/— gM cab ) and (f ab ,N£ b ) with respect to this 
result? At first sight, it might seem that the separation of the S(TS) term into TSS and SST terms 
just corresponds to the separation of the terms with the variation of the metric and the terms with the 
variation of its first derivatives. We have explicitly verified that this is not the case by considering the 
splitting 

S(V=gV c ) = -5{V~gg ab M c ab ) = -M c ab 8f ab - f ab 5M c ab (ill) 

which did not provide us with the TSS + SST splitting. The next natural question would be whether the 
{TSS, SST} structure corresponds to the {pSq, qSp} structure since (g ab , \J~^g~M cab ) and {f ab , \J— gN^ b ) 
are canonically conjugate variables. To answer this question, we looked at another canonically conjugate 
pair (g ab , —y/—gM£ b ) (See Eq. (8)). The corresponding variation 

6(V=gV c ) = -5{g ab ^—gAr ab ) = -^M c ab Sg ab - g ab S(y/^M c ab ) (112) 

also failed to give us the {TSS, SST} splitting, proving that it is not a purely a result of the {pSq, qSp} 
structure. Thus, g ab and f ab and the corresponding canonical momenta are indeed special for our 
purpose. Although we have not yet discovered a completely satisfactory reason as to why this must be 
so, we do have some indication that it must be related to the scaling arguments detailed in Appendix A. 
We hope to return to this issue in a future work. 

From these we can also obtain a thermodynamic identity TdS = dE + dW. Now to obtain this identity, 
we shall borrow some of the results derived in [16]. To begin with, we need to find the variation of \fa. 
If A is the affine parameter corresponding to the tangent vectors of outgoing and ingoing null geodesies, 
then near the horizon we find that 

A = A H + ^Kn 2 + 0(n 3 ) , (113) 

where A = Ah is the location of the horizon. Expressing Eq. (104) in terms of the affine parameter, we 
obtain 

N = V2^(A - A H ) 1/2 + 0((A - Ah?' 2 ); 

<J AB = [a H {v,AH)] AB + K- 1 [a 2 (y,AH)] AB {A-AH) + 0{{A-A H f /2 ) . (114) 
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Let us assume that the transverse metric is independent of whatever parameters are present in the metric. 
Then, the variation of yfa can be written as 

6yfc ee S A ^ = ^-5A = \^o AB ^5A = ±yfc* 2 5A , (115) 
OA 2 oA 2k 

where in the last step Eq. (114) has been used and 02 = o ab \o^\ab- Therefore, we have 

N? k 6f* = ^g 2 SA ; (116) 
Near the horizon, the nn and tt components of the Einstein equations can be written as [16,51] 

\{?2-R\\) =8xTZ =^Tl , (117) 

where R\\ denotes the Ricci scalar on the two-dimensional surfaces of constant n and t. Substituting this 
in Eq. (116), we have 

N? k Sf jk \ H = V^(167tT; + R n )5A . (118) 

Here, * can stand for either n or t because, in our case, T™ = T\. Integrating over the transverse variables 
and using Eq. (110), we obtain 

/ d 2 x^N: b Sf ab = f ^xlv^CC + ^r)5A = TSS. (119) 
167T J J Ion 

The first term in the above equation can be interpreted either as an J PdV term or an J(—p)dV term, 
where P is the transverse pressure and p is the energy density at the horizon, whereas the last term can 
be interpreted as the variation in energy associated with the horizon SE (for details, see [16]): 

&xjLy/a-^8h (120) 

167T 

Then, we arrive at either of the following two equations: 

SE = TSS+ [ pdV; (121) 



SE = TSS - J PdV. (122) 

Thus, we see that the variations of the surface term on a horizon in a general static spacctime can be 
given a thermodynamical interpretation. 

We shall now specialize to a spherically symmetric metric in order to gain some more physical insight 
into our results. 



5.2.2 An application: spherically symmetric metric 

In this subsection, we shall specialize the results of the last section for a spherically symmetric metric of 
the form 

dr 2 

ds 2 = -f(r, \)dt 2 + + r 2 d0 2 + r 2 sin 2 Odcf . (123) 

f{r, A) 

Here, A is a parameter in the system, like the mass M for a Schwarzchild metric. We shall assume that 
there exists a horizon at r = such that f{rh) = 0. In order to use the results in the last section, we 
need to write this metric in the form of Eq. (103). Define a parameter n such that dr? = dr 2 / f(r) and 
n = at the horizon. Near the horizon, we can expand f(r) as f(r) w 2k (r — r/j), where k is the surface 
gravity associated with the horizon at r/j. Taking square root and integrating from r = rh to r = r, 
we obtain n = */2//c(r — r^) 1 / 2 . Comparing with Eq. (113), we see that r plays the role of the affinc 
parameter here. Having obtained this expression, if we now make the identification iV 2 = f(r) and note 
that (Jab is the metric on the surface of the sphere with radius we have made all the connections 
necessary to carry over the results in the previous section. But rather than referring back to the results in 
the last section, we shall rederive these results (in a slightly different manner) since many of the integrals 
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in the last section can be explicitly evaluated in the current case. Also, it will be a good consistency 
check on our results from the last section. 



For the spherically symmetric metric, the radial component of the surface term in the Einstcin-Hilbcrt 
action Aeh (see Eq. (1), Eq. (3) and Eq. (10)), i.e. y/^gV r / 16ir integrated over 9 and <f> at the 2-surface 
at r = r/j, gives us —TS. 

1 



16tt 



d.9 



d4>^gV r = -TS. 



(124) 



We shall now look at the variation of the surface term along the affine parameter, which is the radial 
coordinate in our case, splitting the variation into a qSp part and a pSq part as in the previous section. 
In this section, we shall consider this variation as being due to the variation of the parameter A due to 
which the horizon will undergo a shift in position. The variations 5X and Srh are then connected by the 
relation /(r/j, A) = f(rh + Srh, A + S\) = 0, which implies 



%SX = -fSr h , 

a condition that we shall use to simplify our expressions. The pSq term in this case is given by 



(125) 



1 

16tt 
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2jv 



d<j>N r ab Sf ab 



(126) 



where we have used the relations that Hawking temperature T = k/2h = f /Att and Bckcnstcin-Hawking 
entropy S = A±/4. As can be expected from Eq. (124), Eq. (99) and Eq. (126), we can obtain the qSp 
part as 



1 

167T 



(W 



d(br b sN: b 



Airrl . r ._ d . 3,,/'., 



SST. 



(127) 



Eq. (126) can be rewritten making use of the Einstein equations to provide a clearer physical picture. 
The Einstein equations for the spherically symmetric metric are given by 



Gt s~ir 



Gl = G 



rf + f-l 

4 _ 2f + rf" 
2r 



87rr; ; 



8ttTI = 87rTI 



(128) 
(129) 



Substituting in Eq. (126) from Eq. (128), we obtain 



1 

167T 



dO 



d4>K b 6f 



ah 



^ + T}{Airrl8r h ) = ^ + T r r (^r 2 h 6r h ) = TSS. (130) 



The factor AnrybTh, is just dV, the change in volume when the horizon shifts outward by an amount Sr^. 
Let us assume that the region outside the horizon contains a perfect fluid at rest. The energy-momentum 
tensor will be given by 

Tab = (p + P)u a u b + Pg ab (131) 

and we have Tq = — p and T, r = P. Then, if Srh/2 (which is SM in Schwarzchild case) is interpreted as 
the change in energy SE, Eq. (130) can be written either as 

5E = TSS + p6V (132) 



or 

SE = TSS - P5V. (133) 

Among these two interpretations, Eq. (133) is in the familiar form of the first law of thermodynamics 
but Eq. (132) maybe physically more intuitive as it makes clear that there are two contributions to 
the change in the energy: one contribution from the change in the area of the horizon (TSS term) and 
another from the energy density of the matter engulfed by the horizon when the horizon expands outward 
(pdV term). 
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To understand this interpretation clearly, let us consider the special case of the Reissner-Nordstrom 
metric. We can rewrite Eq. (130) for this case by substituting for Sr^ in terms of SM and SQ using 
r h = r± =M± v/M 2 - Q 2 . We obtain 

,■2 TV 



1 

167T 



de / d<j>N r ab 6f 



lit, 



SM , = (134) 

M ± y/ M 2 - Q 2 



where the left-hand side has already been identified with TSS. Further, with SE = Srh/2, we can 
obtain 

SE - TSS = -0— (Airr h 2 5r h ) (135) 

The energy-momentum tensor for an electromagnetic field, which is acting as the source, is given by the 
expression [54] 

Tv* = ^ (f^F/ - ± g ^F pa F^ (136) 

For the electromagnetic field of a Reissner-Nordstrom metric, there is only one independent non-zero 
component of the field strength tensor, F tr = —Q/r 2 . Evaluating the Tq component at the horizon, we 
see that Eq. (135) can be rewritten as 

5E - TSS = -T^(4nr h 2 Sr h ) (137) 

which leads to SE — TSS = pdV . Thus, we reproduce Eq. (132). Of course, we could also write this 
relation in the form of Eq. (133). But the current version seems easier to relate to as the right-hand 
side in Eq. (135) represents the contribution from the electromagnetic energy density that the horizon 
engulfs as it expands outward. 



5.2.3 Generalization to an Arbitrary Null Surface 



In the context of emergent gravity paradigm, one extensively uses the concept of a local Rindler frame 
and local Rindler horizon [17]. A local Rindler horizon is essentially a patch of the null surface in the 
locally inertial frame. Every local Rindler observer would attribute a temperature and an entropy to the 
local Rindler horizon. Hence it is natural to extend the above analysis for the metric near an arbitrary 
null surface. We shall now discuss this formalism. 



The metric in the neighbourhood of a null surface is given by 

ds 2 = —2radu 2 + 2drdu — 2rj3Adx A du + (iABdx A dx B 



(138) 



where r — corresponds to the null surface. (A detailed construction is presented in [55,56]). As 
usual, we will first calculate all the quantities on an r-constant surface and then take the r = limit. 
The starting point is to find the normal to an r-constant surface, which is given by n a = iVV a r with 
N = (2ar + r 2 /3 2 )~2. The non-vanishing component of n a is n r = N. The surface term on the r = 
constant surface is given by 



A S ur = TF~ i d 3 xVhn r V r , 
16tt J 



where V r = — (l/g)db(gg r ). For the given metric, it turns out that 



V 



1 r 

(J: 



dull + d r {fi(2ra + r 2 p 2 )} + d A (firp A ) 



(139) 



(140) 



Therefore, 



— ^=<9u/i - V^( 2 « + 2rd ra + 2r/3 2 + 2r 2 /3 A d r f3 A ^ 
-dr/j, - y/jLrd A p A - r -^=d A V. ■ 



2ra + r 2 (3 2 



(141) 



20 



Since, the integration variables in Eq. (139) are u and the transverse coordinates x A , it is convenient to 
take the r — limit at this stage. This reduces the above equation to the following form: 



\/hn r V r = d u n — 2J~jla 



Assuming the Taylor series expansion: 



y/Ji = y^WO^) + rf(u,x A ) + 0{r 2 ) 
we have d uy fji = rd u f(u,x A ) + 0(r 2 ), and hence near the null surface we can write: 



Therefore Eq. (142) reduces to \fhn r V r 
find 



dufi = 2d uy /Ji = 0; ^JJi 



/,(°) 



(142) 



(143) 



(144) 



2\J fi^a. Finally, substituting this in Eq. (139), wc 

(145) 



A sur = -— f dud 2 x A \/fi ( -°'> a . 
8tt J 

We next expand a in a Taylor series as 

a(r, u, x A ) — a (x A ) + rg{u, x A ) + . . . . 
Using this in Eq. (145) and then performing the integration from u = to u wc find 

a A ± 



8- 



whcre we defined an average surface gravity «o as 

fd 2 x A ^/JMa (x A ) 



a 



Ax 



(146) 



(147) 



(148) 



Thus the surface term, calculated on the null surface can be interpreted as A sur = —uTS, where 
T = ao/2ir and S = A ± /A. 

Next we will calculate the terms in Eq. (102), which arise due to the metric variation, on the null surface. 
We only need to compute f ah 5N^ b since the the other term can be identified using Eq. (147). Expanding 
f jk SNJ k under the metric (Eq. (138)) we obtain, 



f k SNL = 2f ur SN r ur + f rr 5N; r + 2f rA 8N r rA + f AB 6N r A 



AB 



^d r a + ^p A d r p A + ^n Ac d uf iAc 
,rl 



+ ^JI(2r a + r A f3 2 )S^-fi A ^d rt i AC 



+ 2rp A ^6[~-d r [3 A + -^d r VAC + ^ 1J d A Li 
+ JJl[i ab 8 -d u( i A B + \{P 2 - a)d r HAB ~ D A (3i 



BD 



(149) 



where 



N r 

U 

N' r 

n:. 



_r r + ir a =-(t u - r r + r- 

ur 1 2 au 2 \ uu ur J 

'\d r a + ^p A d r (3 A + \v AC duVAc ; 



-pa 



1 Ar - 2^ °rfJ>AC , 



1 rA + 2 L aA — ^ 



ru j^r I j^B 

uA ~ 1 rA + 1 AB 



N 



AB 



-^d r (3 A + \^p C d r ii AC + ^ BD d A ^BD ; 

' r AB = 7;d u HAB + ^(/3 2 ~ a)d r fi A B - ^(D A 0B + DbPj^J 



(150) 
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with a = — 2ra and Pa = —t~Pa being substituted in the expressions. (To calculate the components 
of Ng c we have used Eq. (16) and the expressions for the connections evaluated in [56].) Near the null 
surface (r = 5r = 0), the above result reduces to 



f jk 5NJ k = 2y/fis[- \d r a+ ^ AC d uf i A c] +VTw AB s[\d u VAB 



(151) 



Finally, use of Eq. (143) and Eq. (146) lead to vanishing of last two terms near r = and thus we find 
p k 5NJ k \ r =o = 2^/J^5a . Therefore, the analogue of Eq. (110) in this case would be 



which can be represented as 



^- J d 3 xf k SNJ k = uSST 



Now since A* 



uTS, the other term in the variation will be 



-!- [ d 3 xNLSf jk = uT5S 
16tt J ]k J 



(152) 
(153) 

(154) 



Thus, we have generalized our thermodynamic interpretation to the case of the variation of the surface 
term in the Einstcin-Hilbert action evaluated over an arbitrary null surface. 



5.3 Connection with the ADM formalism 

In this section, we shall make some comments relating our formalism with the standard ADM formalism 
[57]. We expect to find some parallels between the two, since the ADM formalism also uses the language 
of canonical variables. The major difference, of course, is that we have not assumed a particular foliation 
of the spacetime. In the ADM formalism, we assume a foliation of the spacetime with a family of non- 
intersecting spacelikc hypcrsurfaces. The dynamics is then generally considered in the volume bounded 
by two such spacelike surfaces and two timclike hypcrsurfaces, with one of the timelike hypersurfaces 
assumed to be at spatial infinity. In the following treatment, we shall borrow the required expressions 
from Chapter 12 of [32]. 

The dynamical variables in the ADM formalism are the components of h a p = g a /3, the induced metric 
on the spacelike surfaces of the foliation. The definition of the term "canonical momenta" in the ADM 
context is different from the definition that we have adopted for our formalism in that the canonical 
momenta in ADM formalism refers to the the derivatives of the Lagrangian with respect to the time 
derivatives of the metric components and hence correspond to ^/—gAI aab in our formalism. The ADM 
Lagrangian, Ladm, leads to non-zero canonical momenta conjugate to h a p while the momenta corre- 
sponding to the other metric components vanish. (This does not happen for £ qU ad and one can find 
that M 000 and M 00b are in general non-zero.) The canonical momentum corresponding to h a p is given 

by 

«/3 d 



P" r = a/a , A V^Ladm) = -Vh{K a P - Kh af} ) . (155) 



d{d h aP ) 

Here, K a @ represents the contravariant components of the extrinsic curvature K mn = — /i^V n n . The 
extrinsic curvature for this particular foliation is given by 

K mn = _ Nh am h nb T O ab = N h a ™h nb N% b (156) 

Note that all the indices in the above expression can take only spatial values since h 0m = 0. Taking the 
trace, we obtain 

2VhK = -y^gV - 2^/=gn m n n M 0mn (157) 

which is the expression relating the Gibbons-Hawking- York counterterm [31,34] with the surface term 
of the Einstein-Hilbert action (see e.g. Exe. 6.3 of [32]). We next write the expression for p a @ , given 

by 

= 2V=9-Qt?fT yS = -2^g-Q%fN° s (158) 



(h) Mi-^V-W^) iv 7<5 

■b ^ac^db _ uab^dc ;„ OI ,„l„ mi 

(h) 



where 2Q^ b = h ac h db - h ab h dc , in analogy with Eq. (2) 
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Having thus made the necessary connections between the variables, we can look at the variations of the 
action in the two formalisms. In fact, the variation obtained in Section 12.4.3 of [32] (see Eq. 12.111) is 
similar in structure to the integrated version of Eq. (22) with a qSp surface term. This variation is given 

by 

(159) 



/ &W=gR) = [ rf 4 
Jv J v 



-gG ab Sg ab + / d 3 x eh ab Sp ab 

'V JV JdV 

where e is —1 on spacelikc parts of the boundary dV and +1 on the timclikc parts. Here, h ab is the 
induced metric on the surface of integration and p ab = — y/h(K ab — Kh ab ), where K ab is the extrinsic 
curvature corresponding to that surface. Eq. (155) is a special case in which the surfaces are constant 
time slices. In obtaining this expression from the variation of the action, a surface term has been thrown 
away assuming that the surface of integration is compact. Comparing with the integrated version of 
Eq. (22), we obtain 

- / d 3 x g lk n c 5{^M clk )] = f d 3 x e h ab S P ab . (160) 

JdV JdV 

Here, n c is the unnormalized normal to the integration surface. For example, if we were integrating over 
the upper time slice of the boundary of a usual ADM integration volume, an a; =constant surface, we 
would have n c = —5®, where the minus sign ensures that n c is in the direction of increasing time. The 
normalized normal in this case would be given by n c = —N6®. On the other hand, if our integration 
volume was inside an r = constant surface, we would have n c = 5 r c and the normalized normal would be 
given by n c = (l/g rr )S®. Thus, we have obtained the correspondence for the qSp variation term. 



For connecting up with the pSq variation term, consider the following relation: 



a b 



d 3 x 2VhKe 



i)V 



d A x e p ab 8h 



(161) 



c)V 



If the second term was the variation of our usual surface term in the Einstein-Hilbert action, we could 
have compared the last term directly. with the integral of NSf term or Mdg term. But we have here the 
Gibbons-Hawking- York countcrterm instead of the surface term in Hilbcrt action. To obtain the desired 
relation, we first write down the structure of the usual variation of the Einstein-Hilbert action. This is 
given by 



5(y/=g~R) = [ d 4 x^g~G ab 6g ab + S[ [ d : 

JV \JdV 

Comparing Eq. (162) with Eq. (161), we obtain 



x e n 



\[hV 



d 3 x n t VhM lab Sg ab (162) 



tPx e p ab Sh 



d 3 x mVhM iab 



fig a b + 5 



,)V 



[ d 3 x e VhiriiV* - 2K] 

JdV 



(163) 



We can now use the following result (see Exe. 6.3 in [32]): 

V a n a - 2K = 2h ab d b n a - n m h ns d n g sm = -n m h ns d n g sm (164) 

where we have used the result d b n a = —(e/2)n b n,inj .d a g % ^ and h ab n a = (see Section 12.4.3 in [32]). 
If the metric has no off-diagonal components with respect to the coordinate which labels the surfaces of 
foliation, the right-hand side of Eq. (164) vanishes on the foliation surfaces. In such a case, if we assume 
that the boundaries of our integration volume other than the foliation surfaces do not contribute, (i.e. in 
ADM formalism, for example, assuming the integration region is between two time-slices and a time-like 
surface at spatial infinity where all fields go to zero), we can write Eq. (163) as 



OV 



d 3 x e p ab 6h ab = / d 3 x ni VhM iab 5g ab 



(165) 



This is the desired connection between the p8q variations; but unlike Eq. (160), this relation is valid 
only when the coordinates are chosen such that 2K = riiV 1 , which can be achieved by demanding a 
metric block diagonal with respect to the foliation coordinate. (If we take t = constant surfaces for our 
foliation, for example, then the shift function should vanish.) Further, only the foliation surfaces should 
contribute to the surface integral. Hence, to summarise, we have a relation between the u qSp v variations 
in ADM formalism and our formalism, Eq. (160): 



c)V 



d 6 x g ik n c 6(^M clk )} = d 3 x e h ab Sp 



ab 



dV 



(166) 
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and, in coordinates in which the metric is block diagonal with respect to the foliation coordinate and 
non-zero contributions to the surface term come only from the foliation surfaces, we also have the 
corresponding relation between pSq variations, Eq. (165): 

f d 3 x ni ^/hM iab 5g ab = [ d 3 xe p ab 5h ab . (167) 

JdV JdV 



5.4 Action Principle as a Thermodynamical Extremum Principle 

In the previous analysis, we have shown that £ qua d can be interpreted as a Hamiltonian (Section 4.4) 
- more precisely a Hamiltonian density- and the surface integral arising from C sur leads to — TS on a 
horizon (Section 5.2). This interpretation leads to the interpretation of the Einstein-Hilbert action as a 
free energy, Aeh = tF = t(E — TS), where r is the range of time integration in any static geometry, 
which has been explicitly demonstrated for static metrics in Einstein gravity [24] and also for Lanczos- 
Lovclock models [30] . Wc shall now use this interpretation to formulate the principle of extremisation of 
gravitational action as a thermodynamic extremum principle. 

The variation of the action is given by integrating Eq. (26) over a spacctimc volume. This variation is 
given by 

16nSA EH = [ d i x5{yf^R)= [ d 4 xR ab 8f ab - [ d 3 xVhn i g : > k 8N l jk , (168) 

Jv JV JdV 

where we have rewritten the volume integral of the total divergence as a surface integral. 

Now consider the variation of the action on-shell. For pure gravity (R a b = 0), the above variation reduces 
to 

SA BH = --i- / d 3 xVhn^ k 5N] k = --L / d 3 xf k 6Nf k , (169) 
167rJ av J 167T J 8V ■> 

on the X =constant surface. Here X = n for a static spacetime while X = r for null metric. From 
Section 5.2, we know that the right-hand side, when evaluated on a horizon, can be interpreted as a 
—tSST term. Therefore, using Aeh = t(E — TS) as has been argued, we find that the Eq. (169) can 
be written thermodynamic identity: 

S(E-TS) = -S6T; i.e., SE = T5S (170) 



Next we shall consider the inclusion of the matter action A rn . The usual matter Lagrangians are in- 
dependent of the derivatives of the metric and hence the variation with respect to the metric will not 
involve any surface terms. The variation of the matter Lagrangian on varying the metric is then written 
in the form 

SA m = l f d 4 xV=gT ab 6g ab = ~ f d 4 xV^gTabSf ab . (171) 
* Jv * Jv 

Here, T ab = T ab — (g ab /2)Tl and T ab is the energy momentum tensor corresponding to the matter 
field under consideration. Then, from Eq. (168) and Eq. (171), we can impose the on-shell condition 
Rab = 87rGT a b and obtain 

S[A E H + A m ] = - f d 3 xP k 6N* . (172) 

JdV 

If the matter is perfect fluid, then the matter action, with the on-shell condition, can be expressed 
as 

A m = J v ^d i xP (173) 

where P is the pressure of the fluid [58-60]. For the case of a static spacetime with P independent of f ab , 
the time integration can be performed to give a factor r and the variation of the matter action will reduce 
to tP5V , where V stands for the three-dimensional volume of a time-constant slice. Hence, the left hand 
side of Eq. (172) can be given a thermodynamic interpretation as t8{E — TS) + tPSV. Thus, Eq. (170), 
in a static space time with the inclusion of matter fields in the form of a perfect fluid with P =constant, 
becomes SE = T8S — PSV, which is the thermodynamic identity on the horizon obtained earlier. Wc 
thus see that the variation of the Einstein-Hilbert action allows for a straightforward thermodynamic 
interpretation. 



24 



There is, however, a nicer way of interpreting the gravitational action principle in thermodynamic lan- 
guage, along the following lines. We first note that, if — instead of demanding S[Aeh + A m ] =0 — we 
demand the condition in Eq. (172) we will get the field equations. In such a formulation, we can use any 
spacetime region V and its boundary dV . Consider now a spacetime region bounded by null surfaces. 
Then the surface term on the right hand side of Eq. (172) can be interpreted as giving SST based on 
our earlier result in Eq. (153). Hence, the condition that the surface term vanishes is equivalent to the 
condition that the variations keep the temperature of the null surfaces, as perceived by the local Rindlcr 
observers, constant during the variation. This gives a very direct thermodynamic interpretation of the 
gravitational action principle provided we formulate it in a region bounded by null surfaces. We hope to 
explore this in detail in a future publication. 

6 Conclusions 

The variational principle in general relativity is somewhat peculiar. The key reason is the presence of 
second derivatives of the metric in the Einstein-Hilbert action, which causes the surface term in the 
variation to contain variations of the metric and its derivative. Thus, in order to get the Einstein's 
equation by the usual variational method, we have to fix the metric as well as its normal derivative. The 
main problems with fixing both the metric and the normal derivatives are: (a) The values at which we fix 
the metric and its derivative at the boundaries might not turn out to be compatible with the equations of 
motion derived, (b) If we extend our theory to the quantum domain, we should refrain from fixing both 
the metric and its normal derivative on a spacelike hypersurface to avoid conflict with the uncertainty 
principle. One common method for dealing with this issue is to throw the second derivatives into a 
surface term and then work with the remaining (coordinate dependent) quadratic Lagrangian £ qua d (see 
Eq. (3)), in which case one will have to fix just the variation of the metric on the surface. The other 
common procedure is to use a counterterm to cancel the variation of the surface term [31]. All these 
procedures, one must admit, appear rather contrived. 

In this paper, we discuss an alternate prescription. We may be in better shape conceptually by in- 
terpreting the generally covariant Lagrangian (viz. the Einstein-Hilbert Lagrangian) as a momentum 
space Lagrangian and fixing the canonical momenta on the boundary. We can do this in terms of the 
metric g a b and its canonical momenta with respect to the quadratic Lagrangian, v^M ca& , defined m 
Eq. (7). This suggests that general relativity is better represented as a theory in the space of the canon- 
ical momenta M cab . The immediate direction suggested by this realization is to apply the technique 
of momentum-space path integrals to general relativity, a direction which we intend to explore in the 
future. 

In the process, we discovered the surprising fact that this approach works only with g ab but not with g ab 
and its corresponding canonical momenta. In the case of g ab , the surface variation is found to contain 
the variation of g ab as well as its canonical momenta. We show how this is related to the —d(pq) 
structure of the surface term in the case of the variable g a b and proceed to show how this structure can 
be obtained by simple scaling arguments applicable for homogeneous functions. These arguments also 
allow us to discover another variable f ab = ^/—gg ab , with N£ b representing the corresponding canonical 
momenta. The use of N£ b makes it easy to see that we are actually fixing the variation of the connection 
on the boundary, a fact which needs nontrivial calculation to discover when working with g ab and M cab . 
Further, as already noted in the literature before (sec [36-38,44]), many formulas simplify if we work 
with f ab and N* b . 

The most surprising result of our investigation was the connection between these holographically conju- 
gate variables (HCVs) and thermodynamic quantities pertaining to the null surfaces which act as local 
Rindlcr horizons. The surface term in the Einstein-Hilbert action, when integrated over a horizon, gives 
us the Bekenstein-Hawking entropy of the horizon when the range of time integration is fixed by pe- 
riodicity in Euclidean sector. Without using the Euclidean time integration, and in fact removing the 
time integration altogether, the integral of the surface term over the space variables on the horizon gives 
us [27,28,30] the heat content H sm = TS. In a variation of this integral, it was seen that the term TSS 
is the term with the variation of g a b (or f ab ) and the term with S5T is the term with the variation of 
the corresponding canonical momentum. This result holds: (a)near a horizon in an arbitrary (i.e., not 
necessarily spherically symmetric) static spacetime and (b) near any null surface which acts as a local 
Rindlcr horizon. 
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We believe that this is a very strong result. We know that the variations of T and S are related to 
the variations of the surface gravity k and the area of the horizon respectively. Since n is related to 
the derivative of the goo component of the metric along the direction normal to the horizon, it is clear 
that the variation of the temperature cannot be given by the term with the variation of g ab or f ab and 
must be contained in the term with the variation of the canonical momenta. But it is not clear why 
this canonical momenta term does not contribute to the variation of the entropy To drive this point 
home, we carried out variations with certain other sets of variables (see Eq. (Ill) and Eq. (112) and the 
accompanying discussion) and observed that we do not obtain this separation between the entropy and 
the temperature. Although we do not yet have a clear line of reasoning to offer as to why this must be 
so, we think that the explanation might be related to the scaling properties that were used to arrive at 
the variable f ab . Further, we were able to obtain this result without using the Euclidean time method. 
This seems to suggest that the method of reducing the integral of surface term on a horizon to entropy 
using Euclidean time arguments might not be necessary. 

There are many directions to proceed forward from the work in this paper. It is not clear to us why 
the special variables that we discovered through scaling turn out to be the ones that are related to 
thermodynamic variables. This connection should be further explored. As already mentioned, another 
direction of work would be to try to develop a momentum space path integral approach to general 
relativity. Extension of our results to the case of stationary, or more ambitiously, time-dependent metrics 
and Lanczos-Lovelock models is another obvious line of attack. The simplicity of the scaling argument 
that led us to discover the special nature of the variables g a b and f ab , the naturalness of our prescription 
for the variational approach to general relativity compared to prescriptions existing in the literature as 
well as the intriguing connection with thermodynamic variables on a horizon which is highly unlikely to 
be accidental, suggests that these directions of research would prove to be quite fruitful if pursued. 
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Finding Holographically Conjugate Variables through Scal- 
ing Arguments 



In this section, wc shall try to find alternatives to g a b for which a relation of the form in Eq. (19) holds. 
More explicitly, for Eq. (18) (which we reproduce below): 



q A F A = (A + n)L - di 



1A 



dL 



d{diq A ) 



(A.l) 



we should have q A F = y/—gR and A + \x = 1. We shall restrict ourselves to variables that can be 
obtained by the so-called point transformations [47] i.e transformations where the new variables depend 
on the old variables, but not on their derivatives i.e 



h = h(g ab ,Xi); 



dh 



d{d c g ab ) 



= 0. 



(A.2) 



Further, we shall assume that the transformation has no explicit dependence on spacetime coordinates 
i.e h{g a b,Xi) = h(g a b)- Here we have not specified the index structure of the variable h, but it has to 
be a 2-indexed symmetric object in order to hold the degrees of freedom initially present in the metric. 
First, let us look at the left-hand side of Eq. (A.l). For a point transformation from the variable q to 
the variable s, wc have the result that 



dL 



dt \ dsj ) dsj 



dL 



d_ fdL\ _ dL_ 
dt \ dqi J dsj 



ds.j 



(A.3) 
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Generalizing to our case, we have 



( dL 



dL 



\d{d c g ab ) J dg ab 



( dL \ _ 8L 
\d(d c h) ) ~ ~dh 



dh 
dg ab 



(A.4) 



We shall now assume that h is a homogeneous function of the components of g ab i.e if g ab — > ag ab , then 
h — > a k h for some constant k. If we now contract both sides of Eq. (A.4) with g ab and use Euler's 
theorem for homogeneous functions, we obtain 



9ab 



dL 



dL 



d(d c g ab ) J dg a b 



kh 



0, 



( dL 



\d{d c h) 



dL 
~dh 



(A.5) 



In terms of the notation in Eq. (A.l), we can state this result as follows. If we transform from a set 
of variables qA to another set j a , such that the j 'a do not dependent on the derivatives of qA and are 
homogeneous of degree k in qA, then 

q A F q A — > kf A Ff, (A.6) 

where we have used the subscript in the Euler-Lagrange function to denote the variable that has been 
used. 

Next we shall look at the right-hand side of Eq. (A.l). Under a point transformations, it is easy to see 
that the value of A, the degree of derivatives of the variables, remains a constant. A fallacious argument 
for finding the change in (i is the following, j a being homogeneous in qA of degree k and L being 
homogeneous in qA of degree /j,, a change qA — > aqA would correspond to a change Ja — > erf a and a 
change L — > a^L. Therefore, under a change Ja — > oiJa, we should have qA — > aX'^qA and L — > a^l k L 
and we can conclude that the degree of L in the variable Ja is [ijk. But as is clear from Section 3.2, this 
argument fails for the transformation g ab — > g ab . To find out the error in the above argument and to derive 
the correct result, consider a general term in the Lagrangian of the form (qA)^ {diqs) X ■ The notation 
(<?a) m here corresponds to a term of the form {(qAqs----) H terms} and {diqB) x corresponds to a term 
of the form {(diqAdjqs----) — > A terms}. In terms of the new variables /a, this term becomes 



{q A y{d iqB f = {q A [f c ]Y(^f- di f c ) X , 

dfc 



(A.7) 



where we have made use of the fact that qA does not depend on the derivatives of fs- Now, this 
assumption also tells us that (dqB/dfc) is a function of Ja alone and not its derivatives. Thus, the 
factor containing derivatives of Ja is (dfc) X confirming that the degree in derivatives does not change 
under the transformation. The factor that depends only on the variables j a and not on its derivatives 
is 



{q A [fc]Y{^) X . 
ojc 



(A. 



We have considered Ja to be a homogeneous function of qs, and qs only, of degree k. Since we are 
conserving the degrees of freedom, we should be able to invert these functions and express qs in terms 
of /a as homogeneous functions of degree 1/k in /a- Then, (dqB)/(dfc) will be a homogeneous function 
of Ja of degree {1/k) — 1. Hence, the above function will be a homogeneous function of /a alone, and of 
degree + [(1/fc) - 1]A. 

It is straightforward to generalize our results for (qA)^ '{dqB) X to L and conclude that, under a trans- 
formation from variables qA to variables fs homogeneous of degree k in qA and independent of the 
derivatives of qA, we shall have 



A; A^A 



A 



/' 



X + fJt 



(A.9) 



Thus, from Eq. (A.l), Eq. (A.6) and Eq. (A.9), the form of Eq. (19) is conserved only if we transform 
to a variable which is homogeneous of degree k = 1 in g ab . This obviously breaks down for g ab , which 
has k = — 1. The determinant of g ab , g, does not have the required number of degrees of freedom and 
has k = 4 in four dimensions and is also unsuitable. But the variable 



f ab = V~gg ab 

has the required number of degrees of freedom and has k 
contravariant in its two indices, to g ab . 



1, thus providing a useful alternative, 



27 



B Useful properties and relations pertaining to / 



ah 



In this appendix, wc shall list out some useful properties and relations pertaining to f ab = \/—gg ab - Its 
determinant is given by 

./ = dct(/ ab ) = dct(g ab ) = g. (B.l) 

Hence, all the s/^g factors that hang around in expressions can be replaced by y/—f. In order to use 
the well-known formulas used in variation of the gravitational action while working with f ab , we need to 
relate the variation of g ab to the variation of f ab . This relation is given by 



Sg 



lm 



<\f 



J Jab _ £> ab <JJ 



-/ 



(B.2) 



where B l £ = (l/2)(6 l a 5^ + 8[5™) - {l/2)g lm g ab = 5[ a 6™ - (l/2)g lm g ab . We shall take B l ™ to be B lm ab 
so that Bi ma b = gi a g m b - {l/2)gi m g a b- B 1 ^ satisfies the relation 



T>lm r>ab cZ rm 
U ab B cd = d (c°d)- 



(B.3) 



This relation is valid even if we remove the explicit symmctrisation in both B l b m and the right-hand side of 
this relation. Using this relation, we can easily invert Eq. (B.2) and obtain (see Chapter 6 in [32]) 

Sf lm = V=gB l a " b l Sg ab . (B.4) 

Hence, we have, for any two-indexed object X\ m 

X lm Sg lm = -^[X ab -\g ab g lm X lm ]5f ab 

V^g 2 

= -L[X a6 - \g ab X].5f ah (B.5) 

Therefore, 



V=gG ab 5g ab = R ab Sf ab . 
The Euler-Lagrange function (see Section 4.3) F ab for the variable f ab is 

d(y/^gL qulid ) q r 9(^gL quad ) 1 



F„ 



Rab 



df ab CL d(d c f ab ) 

This expression can be arrived at either by explicit computation or by staring at Eq. (26). 
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C Proof of the Conservation Equation di{t l k — 167T yf^gTfy = 

Consider the object &it\. We have, using the definition in Eq. (90), 
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= d t 
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- U 
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d k f ab 



(C.l) 



It is clear from the above equation that, in the absence of matter, the object t\ will be conserved once 
we impose the equations of motion. 

To generalize to the case with matter present, start from the Bianchi identity ^/— gSI t G\. = 0. Then, we 
have ^f r gV l R l k = (^g/2)V k R, implying 



2d i (V=?fl£) = V=99 ab VkRab + V=gR ab d k g a b = f ab d k R ab 



(C.2) 



2<S 



(C.3) 

Thus, we obtain the general 
(C.4) 
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